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ABSTRACT 

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are pre- 
sented for the first time and analyzed using a rate equation model. The energy barriers for the relevant 
diffusion and desorption processes are obtained. They turn out to be significantly higher than those 
obtained earlier for polycrystalline silicates, demonstrating the importance of grain morphology. Us- 
ing these barriers we evaluate the efficiency of molecular hydrogen formation on amorphous silicate 
grains under interstellar conditions. It is found that unlike polycrystalline silicates, amorphous silicate 
grains are efficient catalysts of H2 formation within a temperature range which is relevant to diffuse 
interstellar clouds. The results also indicate that the hydrogen molecules are thermalized with the 
surface and desorb with low kinetic energy. Thus, they are unlikely to occupy highly excited states. 

Subject headings: dust — ISM: molecules — molecular processes 
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1. INTRODUCTION 

H2 is the most abundant molecule in the interstellar 
medium (ISM). It plays a crucial role in the initial cool- 
ing of clouds during gravitational collapse and is involved 
in most re action schemes that produce other molecules 
(jDulev fc W illiams 1984). It is widely accepted that 
H2 formation in the IS M takes place on dust grains 
(jCould fc Salpeter 1963f ). In this process, H atoms that 
collide with a grain and adsorb on its surface quickly 
equilibrate and diffuse either by thermal activation or 
tunneling. They ma y encounter each other and form H 2 
molecules (iWilhams 1968t iHollenbach fc Salpeter 1970f ). 
or desorb thermally in atomic form. 

In recent years, we have conducted a series of ex- 
periments on the formation of molecular hydrogen 
on du st grain analogues suc h as polycrystalline sil i- 
cates (jPirronello et al. 1997al: iPirronello et al. 1997bl ). 
amorphous carb on (iPirronello et al. 1999|) and amor- 
phous water ice (|Mamco et al. 20011: iRoser et al. 2002t 



Dulieu et al 20051 



Amiaud et al. 2006t IWilliams et al. 



iRoser et al. 2003t iPerets et al. 20051 ). under astro- 



physically relevant conditions. In these experiments, 
the surface was irradiated by beams of H and D 
atoms. The production of HD molecules was mea- 
sured during the irradiation and during a subsequent 
temperature programmed desorption (TPD) experi- 
ment. To disentangle the process of diffusion from 
the one of desorption, separate experiments were 
carried out in which molecular species were irradiated 
on the sample and were later induced to desorb. 
Related studies were done on amorphous ice sur- 
faces ()Hornekaer et al. 20031 iHornekaer et al. 2005t 
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20071) . 

The results were analyzed using rate equation models 
and the energy barriers for the relev ant diffusion and 
desorption processes were obtained (iKatz et al. 19991 : 
ICazaux fc Tielens 2004 IPerets et al. 20051 ). Using these 
parameters, the conditions for efficient H2 formation on 
different astrophysically relevant surfaces were found. In 
particular, the formation of H2 on polycrystalline sili- 
cates was found to be efficient only in a narrow temper- 
ature window below lOK. Since the typical dust grain 
temperature in diffuse interstellar clouds is higher than 
lOK, these results indicated that polycrystalline silicate 
grains cannot be efficient catalysts for H2 formation in 
most diffuse clouds. 

In this Letter we present, for the first time, experi- 
ments on molecular hydrogen formation on amorphous 
silicates and an alyze the results u sing a suitable rate 
equation model (jPerets et al. 20051 ). Using the param- 
eters that best fit the experimental results, the effi- 
ciency of hydrogen recombination on grains is obtained 
for a range of conditions pertinent to diffuse interstellar 
clouds. It is found that unlike the polycrystalline silicate 
grains, amorphous silicate grains, which are the mai n 
silicate component in interstellar clouds ijTielens 20051) . 
are efficient catalysts for H2 formation within a broad 
temperature window that extends at least up to about 
14K. 

2. EXPERIMENTAL METHODS AND RESULTS 

The apparatus consists of an ultra-high vacuum cham- 
ber housing the sample holder and a detector. The sam- 
ple can be cooled by liqiud helium to ~5K, as measured 
by a calibrated silicon diode and thermocouple placed in 
the back of the sample. A heater in the back of the sam- 
ple is used to maintain a set temperature between 5 and 
30K during the irradiation phase of the experiment. The 
sample and detector can rotate around the vertical axis. 
Prior to a series of measurements, the sample is heated 
to 380-400K. During a series of measurements, the sam- 
ple is taken periodically to 200-250K to desorb condens- 
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ables. Hydrogen and deuterium gases are dissociated in 
two radio-frequency dissociation sources, with dissocia- 
tion effi-ciency of 80-90%, and are sent into the sample 
chamber via two tr iple differentially pumped beam lines 
(jVidali et al. 20 041 

In the experiments reported here, we used beams of 
low fluxes and short dosing times. Using the standard 
Langmuir-Hinshelwood analysis , plotting the total yield 
of HD vs. the exposure time (i Biham et al. 200ll ). we 
estimated the coverage to be a small fraction (a few per- 
cent) of a monolayer (ML). This is still far from inter- 
stellar values but is within the regime in which results 
can be safely extrapolated to diffus e cloud conditions 
(|Katz et al. 19991 : iPerets et al. 200^ . The interstellar 
dust analogues we used are amorphous silicate samples, 
(Fe^;, Mgi_^)2Si04, prepared by one of us (J.R.B.) by 
laser ablation (wavelength 266 nm) of a mixed MgO, FeO 
and Si02 target in an oxygen atmosphere (10 mbar). 
The optical and stochiometric characterization of sam- 
ples produced with this technique is given elsewhere 
(|Brucato et al. 20021 ). The results reported here are for 
a sample with x = 0.5. 

The experiment consists of adsorbing hydrogen atoms 
onto the surface while monitoring the amount of hydro- 
gen molecules that are formed. To increase the signal to 
noise ratio, hydrogen and deuterium atoms are used and 
the formation of HD is monitored. The measurement 
of HD formation is done in two steps. First, we record 
the amount of HD that forms and comes off the surface 
while the sample is being dosed with H and D atoms 
(the irradiation phase). Next, after dosing is completed, 
in a TPD experiment, the surface temperature is raised 
rapidly and the rate of HD desorption is measured (the 
TPD phase). By far, the main contribution comes from 
the TPD phase. 

Irradiations with beams of H and D ("H+D" there- 
after) were done on an amorphous silicate surface in or- 
der to explore the formation processes of HD molecules. 
The H-l-D irradiation runs were performed with differ- 
ent irradiation times (15, 30, 60, 120 and 240 s), at a 
surface temperature of Tq ~ 5.6K. In a separate set of 
experiments, beams of HD molecules were irradiated on 
the same surface. During the TPD runs, the sample 
temperature was monitored as a function of time. The 
temperature ramps T{t) deviate from linearity but are 
highly reproducible (see inset in Fig. 2). 

The desorption rates of HD molecules vs. surface tem- 
perature during the TPD runs are shown in Fig. 1, 
for H-|-D irradiation on polycrystalline silicate {circles) 
and amorphous silicate (squares) surfaces, with irradi- 
ation times of 120 s. The TPD curve following irra- 
diation of HD molecules on an amorphous silicate sur- 
face is also shown (crosses). The results of current ex- 
periments of H+D irradiation on amorphous silicates, 
clearly differ from those of earlier experiments on poly- 
crystalline silicates. The desorption curves from amor- 
phous silicates contain two wide peaks, located at a 
significantly higher temperatures than the single nar- 
row peak obtained for the polycrystalline silicate. The 
higher peak temperatures indicate that the relevant en- 
ergy barriers are larger, while their large width reflects 
a broader distribution of the energy barriers of the HD 
desorption sites. The TPD curve of HD desorption from 
amorphous silicates, after irradiation with HD molecules 
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Fig. 1. — TPD curves of HD desorption after irradiation of 
HD molecules (x) and H-l-D atoms (□) on amorphous silicate. 
Also shown, for comparison, is HD desorption after irradiation with 
H-l-D atoms on polycrystalline silicate (o). 
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Fig. 2. — TPD curves of HD desorption after irradiation with 
H-l-D atoms on amorphous silicate with irradiation times of 15 (*), 
30 (x), 60 {□) 120 (-I-) and 240 (o) s. The solid lines are fits 
obtained using the rate equation model. The temperature ramps 
(shown in the inset) are identical for all the runs. 



(crosses in Fig. 1), is qualitatively similar to the curve 
obtained for H-fD irradiation. In particular, the peak 
temperatures are the same. The relative weights of the 
high temperature peaks vs. the low temperature peaks 
are somewhat different. Also, in similar experiments 
with higher values of Tq, a third peak was observed at 
higher temperatures. We attribute this behavior to dif- 
fusion of HD molecules, which grad ually migrate from 
shallow into deep adsorption sites (iPerets et al. 20051 : 
iDulieu et al 20051 : lAmiaud et al. 26m . 

In Fig. 2 we present a series of TPD curves for HD des- 
orption after irradiation with H+D atoms on an amor- 
phous silicate surface with different irradiation times. 
Each curve exhibits a large peak at a lower tempera- 
ture and a broader peak (or a shoulder) at a higher tem- 
perature. The position of the high temperature peak is 
found to be independent of the irradiation time, indicat- 
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ing that this peak exhibits first order kinetics. The low 
temperature peak shifts to the right as the irradiation 
time is reduced. Since in the low-coverage regime stud- 
ied here, the activation energies and the pre-exponential 
factor are not expected to depend on the coverage, these 
results indicate that the low temperature peak exhibits 
second order kinetics. This is unlike the case of irradi- 
ation with HD molecules, where the lower-temperature 
peak does not shift, showing first order kinetics. 

3. ANALYSIS OF THE EXPERIMENTAL RESULTS 

In the model used here, there is no distinction between 
the H and D atoms, namely the same diffusion and des- 
orption barriers are used for both isotopes. Hydrogen 
atoms that stick to the surface hop as random walkers 
and may either encounter each other and form molecules, 
or desorb from the surface. As the sample temperature 
is raised, both the diffusion and desorption rates quickly 
increase. If a large fraction of the energy released when 
two H atoms recombine is transformed into kinetic en- 
ergy of the formed molecule, it would immediately des- 
orb from the grain surface in a high ro-vibrational state, 
and with large tran slational energy. H owever, both our 
experiments on ice ()Perets et al. 2005[) and the current 
experiments indicate that such prompt desorption does 
not occur on amorphous surfaces [but see Perry and Price 
(2003), Tine et al. (2003) and Creighan et al. (2006) for 
other, more ordered, surfaces]. Instead, the newly formed 
molecules dissipate their energy, probably through mul- 
tiple collisions with the rough surface or internal pores. 
These molecules thermalize with the surface and become 
trapped in adsorption sites before they thermally desorb. 
Consequently, the desorbed molecules are not highly ex- 
cited, and desorb only with a thermal energy comparable 
with the grain surface temperature. Therefore, in our 
model it is assumed that the newly formed molecules do 
not promptly desorb, but are trapped in adsorption sites 
with a range of potential barriers. 

The experimental results were fitted using the rate 
equation model described in Perets et al. (2005). The 
parameters for the diffusion and desorption of hydrogen 
atoms and molecules on the amorphous silicate surface 
were obtained. These include the energy barrier for 
the diffusion of H atoms and the barrier for their 
desorption. The value obtained for the desorption bar- 
rier should be considered only as a lower bound, because 
the TPD results are insensitive to variations in E^^, as 
long as it is higher than the reported value. The desorp- 
tion barriers of HD molecules adsorbed in shallow (lower 
temperature peak) and deep (higher temperature peak) 
sites, are given by E^_^{j), where j = 1, 2, respectively. 

The rate equation model is integrated using a Runge- 
Kutta stepper. For any given choice of the parameters, 
one obtains a set of TPD curves for the different irradia- 
tion times used in the experiments. The actual tempera- 
ture curve of the sample, recorded during the experiment, 
is used in the analysis (see inset in Fig. 2). In the first 
step, the barriers E^^{j), j = 1, 2, for the desorption 
of molecules are obtained using the results of the ex- 
periments in which HD molecules are irradiated on the 
surface. To obtain better fits, we incorporate suitable 
Gaussian distributions of energy barriers around these 
two values. In the second step, the barriers for diffusion 
and desorption of H atoms are obtained, using the model 
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Energy barriers for diffusion and desorption. 
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Polycrystalline Silicate 


25 




32 27 


Amorphous Silicate 


35 




44 35, 53 



to fit the results of H-l-D irradiation experiments (Table 
!)• 

The second order behavior of the low temperature peak 
in the H-|-D irradiation experiments can be explained as 
follows. Most HD molecules are formed only when the 
surface temperature is sufficiently high to enable signifi- 
cant mobility of H and D atoms. At this temperature, the 
shallow adsorption sites cannot retain the newly formed 
molecules adsorbed in these sites, which quickly desorb 
from the surface by thermal activation. However, those 
newly formed molecules which are trapped in deeper sites 
do not desorb yet, and remain on the surface until its 
temperature increases further. Thus, the high tempera- 
ture peak exhibits first order kinetics and is located at 
the same temperatures in both HD and H+D irradiation 
experiments. Note that we cannot exclude the possibil- 
ity that some fraction of the molecules desorbed in the 
high temperature peak are formed during irradiation and 
remain trapped in deep adsorption sites until the tem- 
perature becomes sufficiently high for them to desorb. 

At longer H-l-D irradiation times, more atoms are ad- 
sorbed on the surface, and they find each other more 
easily. The low temperature peak (corresponding to HD 
desorption from shallower sites) shifts to lower tempera- 
tures with longer irradation times, thus showing second 
order behavior. This behavior must saturate with even 
longer irradiation times, when the temperature needed 
for HD formation becomes lower than the temperature 
needed for thermal desorption from the shallower HD 
adsorption sites. The HD molecules then form at low 
temperatures, but are trapped in the shallow sites. They 
then await until the temperature is further increased, 
much like the adsorbed HD molecules in the HD irradi- 
ation experiments. Therfore, at longer irradiation times 
the low temperature peak should saturate into first-order 
like behavior and the TPD curves of the H+D and HD ir- 
radiation experiments should become more similar. This 
behavior is confirmed by the TPD curves for H-l-D and 
HD irradiation shown in Fig. 1. 

Using the parameters obtained from the experiments 
we now calculate the recombination efficiency of hydro- 
gen on amorphous silicate surfaces under interstellar con- 
ditions. The recombination efhciency is defined as the 
fraction of hydrogen atoms adsorbed on the surface which 
come out as molecules. In Fig. 3 we present the recom- 
bination efficiency vs. surface temperature for the amor- 
phous silicate sample under flux of 5.2 x 10"^" (ML s~^). 
This flux is within the typical range for diffuse interstel- 
lar clouds, where bare amorphous silicate grains are ex- 
pected to play a crucial role in H2 formation. This flux 
corresponds to gas density of 10 (atoms cm~^), gas tem- 
perature of lOOK and density of 7 x 10"'^^ adsorption sites 
per cm^ on the surface, obtained using the procedure 
described in Biham et al. (2001). 

A window of high recombination efhciency is found 
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Fig. 3. — Calculated recombination efficiency of Hydrogen at 
steady state on amorphous silicate (solid line) and polycrystalline 
silicate (dashed line) vs. temperature, using the parameters ob- 
tained from the TPD experiments. 



between 8-13K, compared to 6-lOK for polycrystalline 
silicate under similar conditions. For gas density of 
100 (atoms cm~^), the high efficiency window for 
the amorphous silicates surface shifts to 9-14K. At 
higher temperatures atoms desorb from the surface be- 
fore they have sufficient time to encounter each other. 
At lower temperatures diffusion is suppressed and the 
Langmuir-Hinshelwood mechanism is no longer effi- 
cient. Saturation of the surface with immobile H atoms 
might render the Eley-Rideal mechani sm more efficient 
in producing som e recombination (jKatz et al. 19991 : 
iPerets et al. 20051 ). Our results thus indicate that re- 
combination efficiency of hydrogen on amorphous sili- 
cates is high in this temperature range, which is rele- 
vant to interstellar clouds. Therefore, amorphous sili- 
cates seem to be good candidates for interstellar grain 
components on which hydrogen recombines with high ef- 
ficiency. 

4. DISCUSSION AND SUMMARY 

The analysis of the TPD curves from amor- 
phous silicate surfaces shows that the relevant 
energy barriers on these surfaces are signifi- 
cantly higher t han on polycrystallin e silicates 
(|Katz et al. 19991 : ICazaux fc Tielens 200l) . A sim- 
ilar tren d was observed in am orphous and porous ice 
surfaces (jWilhams et al. 2007D . These results confirm 
the effect of surface morphology on the distribution 
of energy barriers. This effect can be parameterized 
using a model that provides a quantitative connec- 
tion between the rough ness and the energy barriers 
(jCuppen fc Herbst 20051 ). Our results are consistent 
with this model and can be used to quantitatively con- 
strain its parameters. More specifically, we find a 1.4-1.5 
times increase in the energy barriers of the amorphous 
silicate vs. the polycrystalline silicate surface. This 
gives rise to shifting and broadening of the temperature 
window in which H2 formation is efficient, by a similar 
factor. Scanning electron microscope images show that 
the morphology of the amorphous silicate samples is very 
rough (H. Cuppen, private communication). These sam- 
ples are made of a broad (log-normal like) distribution 



of spheres and agglomerates, from several nanometers 
up to a micron size. It is safe to conclude that our 
present results, together wi t h those on amo rphous 
carbon (jPirronello et al. 19991 : iKatz et al. 1999f ), show 
that amorphous silicate and carbon grains are efhcient 
catalysts for the formation of molecular hydrogen in 
diffuse clouds. However, the results indicate that surface 
roughness is unlikely to extend the window of efficient 
recombination to temperatures of the order of 30-50K 
observed in photo-dissociation regions (PDRs). We 
thus conclude that surface roughness, by itself, does not 
explain the high abundance of H2 in PDRs. 

In the model used here, it is assumed that II2 molecules 
do not desorb immediately upon formation. Instead, 
they stay trapped in the adsorption sites or hop be- 
tween them until thermal desorption takes place. Con- 
sequently, one needs to consider mechanisms for the 
dissipation of the excess energy acquired from the re- 
combination process in order to prevent prompt des- 
orption. For amorphous, porous ice surfaces it was 
shown, using time of flight (TOF) measurements, that 
the kinetic energy of the desorbed molecules is small 
(~3 meV), n amely the excess energy is absorbed by 
the surface (iRoser et al. 20031 : iHornekaer et al. 20031 : 
iHornekaer et al. 20051) . Although we do not have direct 
measurements of the TOF of the HD molecules desor- 
ped from the amorphous silicate surface, the similarity 
between the TPD curves obtained after HD and H-l-D 
irradiations indicates that newly formed HD molecules 
reside and then desorb from the same adsorption sites as 
HD molecules irradiated on the surface. In light of the 
results on both ice and amorphous silicates, it is likely 
than H2 molecules formed on realistic interstellar dust 
would have low kinetic energy and would probably not 
occupy excited vibrational or rotational states. 

In summary, we have analyzed a set of TPD experi- 
ments on molecular hydrogen formation and desorption 
from amorphous silicate surfaces under conditions rele- 
vant to interstellar clouds. Fitting the TPD curves by 
rate equation models, the essential parameters of H2 for- 
mation on amorphous silicate surfaces were obtained. 
These parameters include the energy barrier for diffu- 
sion of H atoms as well as their barrier for desorption 
(considered as a lower bound) . The distribution of bar- 
riers for desorption of H2 molecules is also obtained. In- 
terestingly, a single type of adsorption site for hydrogen 
atoms is identified, vs. two types of sites for molecules. 
The fraction of the adsorption sites, which belong to 
each of the two types is also evaluated. The rate equa- 
tion model provides a unified description of several first 
and second order processes. It enables us to extrapo- 
late the production rate of H2 molecules from labora- 
tory conditions to astrophysical conditions. It thus pro- 
vides a quantitative evaluation of the efficiency of amor- 
phous silicate surfaces as catalysts in the formation of 
H2 molecules in interstellar clouds. We find that the re- 
combination efficiency strongly depends on the surface 
temperature. In particular, the amorphous silicate sam- 
ple studied here exhibits high efficiency within a range 
of surface temperatures which is relevant to diffuse in- 
terstellar clouds. The comparison of the current results 
with earlier ones on polycrystalline silicate surfaces shows 
the importance of surface morphology in molecular hy- 
drogen formation. The results are in agreement with 
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theor etical predictions on th e effects of surface rough- 
ness (|Cuppen fc Herbst 2005D . The results also indicate 
that on amorphous surfaces, newly formed H2 molecules 
are thermalized on the surface and do not promptly des- 
orb. Consequently, H2 molecules formed on and desorbed 
from realistic amorphous interstellar dust are expected 
to have very low kinetic energy and would probably not 



occupy excited vibrational or rotational states. 
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